February 7, 2001 



Anyons and the elliptic Calogero-Sutherland 

model 



Edwin Langmann 



Abstract 

We obtain a second quantization of the elliptic Calogero-Sutherland 
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^jQ ' (eCS) model by constructing a quantum field theory model of anyons 

on a circle and at a finite temperature. This yields a remarkable iden- 
tity involving anyon correlation functions and providing an algorithm 
for solving of the eCS model. The eigenfunctions obtained define an 
elliptic generalization of the Jack polynomials. 

^ ' There is an interesting class of integrable many body systems in one 

dimensions which is known as Calogero-Moser-Sutherland systems JTJ, [2|, [3|, |J. 
*5 ■ These systems describe an arbitrary number of identical particles interacting 

£H ! with a two-body potential which, in the general CcLSG, IS cL Weierstrass elliptic 

function. Apart from their purely mathematical significance, these systems 
.^J ; are also of considerable physical interest as they are relevant to remarkably 

^ ■ many different topics including fractional statistics and anyons, quantum 

chaos, QCD, and two dimensional quantum gravity (reviews of the first two 
resp. three latter topics are, e.g., B resp. [0]). 

In the important limiting cases where the elliptic two-body potential be- 
come rational and trigonometric, algorithms to solve the quantum version of 
these systems were discovered by Calogero QTJ and Sutherland 0, and the 
mathematical properties of the solutions thus obtained have been studied ex- 
tensively (see e.g. J?], |J). In this letter we give a concise description of various 
novel results on the elliptic Calogero-Sutherland (eCS) model, culminating 
in an algorithm which provides a generalization of the above-mentioned so- 
lution of Sutherland (detailed proofs with technical details will be published 
elsewhere || |H]]). In particular, we find a second quantization of the eCS 
model as quantum field theory model of anyons on a circle and at finite 
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temperature, generalizing recent results on the Sutherland model and zero 
temperature anyons [[11], [L2| . We then use this second quantization to prove 
a remarkable identity [Eq. ([19]) below, which is a particular regularization of 
the identity in Eqs. (|34])-(|35|)1 which finally yields an algorithm to construct 
eigenfunctions and eigenvalues of the eCS model. In the trigonometric limit 
our algorithm simplifies to one which is equivalent to Sutherland's (i.e., it 



is different but equally simple and yields the same solution) [|10j. To our 
knowledge, solutions of the eigenvalue equation of the eCS Hamiltonian were 
previously known only for integer values of the parameter A defined in Eqs. 



([]]) and ([13|) below ]H| [b|, [L5| , whereas we do not impose such a restriction. 

It is worth noting that in the two particle case, the eigenvalue equation of 
the eCS is also known as Lame's equation which was studied extensively at 



the end of the 19 century (see |]13|) and more recently in [[T4J]. Results on the 



iV-particle generalization of this where obtained in |TjJ using quantum field 
theory techniques which, however, are different from ours. Our approach was 
inspired by the relation of the Sutherland model and the fractional quantum 
Hall effect (FQHE) |16|]: it is known that the edge excitation in FQHE states 
can be described by a ID quantum field theory model of composite fermions 
T7| which is identical with the anyon model used to construct the second 



quantization of the Sutherland model |TT], |12| at odd-integer values of the 
parameter A where the anyons become (composite) fermions [see Eq. (||) 
below] . The second quantization mentioned above [JlTJ, [L2| therefore provides 
a direct link between the Sutherland model and the FQHE. In this letter we 
find that the generalization of the former to the elliptic case corresponds to 
a natural generalization of the latter: going from zero- to finite temperature. 
We now fix our notation. The eCS Hamiltonian is defined by 





j=l 3 j<k 

where —it < Xj < tt are coordinates on a circle, 7 > —1/2 is the coupling 
constant, j,k — 1, . . . , N, and the interaction potential is 

V ( r ) =-^2 lo S^(H ( 2 ) 

with the Jacobi theta function -$i(r) = const, sin(r) n^°=i[l — 2g 2n cos(2r) + 
q 4 ] flI8| . Note that V(r) is equal, up to an additive constant, to the Weier- 
strass elliptic function p(r) with periods 2ir and i/3 where q = exp(— (3/2) 
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The particle number N is arbitrary and (3 > 0. The Sutherland model 
corresponds to the limiting case q = where V(r) = (1/4) sin~ 2 (r/2). Of 
course, the differential operator in Eq. ([[]) does not specify a unique quantum 
mechanical Hamiltonian (i.e., it has many different self-adjoint extensions), 
but our approach will automatically specify a unique self-adjoint operator 
which, for q — 0, coincides (essentially) with the one diagonalized by Suther- 
land fli~0| . (This self-adjoint extension is particularly 'nice' and, for 7 > 0, 
essentially the Friedrich's extension |19| . 



We now describe the construction of anyons using ID chiral bosons [17] 
which can be made mathematically precise using the representation theory 
of the loop group Mapi^S 1 ; U(l)) [Q. We recall that formally, an anyon field 
4>(x), parametrized by a coordinate — it < x < tt on the circle, is characterized 
by the exchange relations (p(x)(f>(y) = e ±mX (f)(y)(f>(x) where A > is the so- 
called statistics parameter. To construct such anyons we start with boson 
operators p(n), n integers, together with an invertible operator R and obeying 
the relations 

[p(m), p(n)) = m6 m „ n , [p(n), R] = 5 nfi R . (3) 

It is natural to interpret Q = p(0) as charge operator and R as charge 
rising operator. The standard representation of this algebra is an irreducible 
highest weight representation on a fermion Fock space (see e.g. [|l2j). The 
representation we use is different and will be specified further below. In 
whatever representation, one can define operators 

KJx) = V —pin) e [nx e" |n|£ , (4) 
in 

and then the anyons field can be defined as 

foix) = e" iAQ ^ 2 R e~ iXQx ' 2 I e -^ K ^ 2 (5) 

where e > is a regularization parameter and x • * means normal ordering 
which depends on the representation and will be specified below; for now it is 
enough to know that x • x amounts to multiplication with some e-dependent 
constant. We stress that introducing this parameter e is a convenient tech- 
nical tool for us which takes care of all the (ultraviolet) divergences which 
otherwise would appear: for e > the quantum field <f> e (x) are well-defined 
operators which can be multiplied without difficulty. Eventually we are in- 
terested in the limit e j which is singular (since the Mx) = <fio(x) are 
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operator valued distributions), but we will be able to take this limit at a 
latter point without difficulty. Using Eq. (j^) and the Hausdorff formula (i.e., 
e a e b = e~f a '^ e b e a for operators a and b such that [a, b] commutes with a and 
b), we obtain exchange relations 

&(x)Mv) = e-™ Xs ^>^My)M*) (6) 

where sgn £ (x) = (1/tt)[x + ^ n ^ e mx ~^ e / (in)] is a regularized sign function 
on the circle. This shows that the <p e (x) is indeed a regularized anyon field 
operator with statistics parameter A. 

We now specify the representation we are using and construct our rep- 
resentation space T . For that we introduce two commuting copies of the 
algebra in Eq. (|), i.e., [/^(m), Pfi(ri)] = m5 m - n 5A,B for A, B = 1,2, and 
similarly for -Ri,2- Then T is the Hilbert space generated by these opera- 
tors from a highest weight vector (vacuum) Q such that p^(n)fi = and 
PA{n)* = f>A(—n) for all n > and that the Ra are unitary operators such 
that < fl, R™VL >= 5 m fl for all integers m (* is the Hilbert space adjoint and 
< -, • > the inner product). We then set 

p(n) = c\ n \pi(n) + s\ n \p 2 {-n) Wn ^ 0, (7) 

and p(0) = pi(0) and R = R\, and it is easy to see that this gives a represen- 
tation of the relations in Eq. (S) provided that (? n — s\ — 1 for all n — 1, 2, — 
In particular we choose 

with q = exp (—(3/2) and /3 > 0. One can prove that this is (essentially) the 
finite temperature representation for the many particle Hamiltonian Tio = 
* P(.~ n )p(, n )/2 x constructed by the usual trick of doubling the degrees 
of freedom (see e.g. |p0|), and the parameter (3 corresponds to the inverse 
temperature |9j . To define the normal ordering prescription for anyon opera- 
tors we determine the creation- and annihilation parts of the operator K e (x) 
such that K = K+ + K~, K~VL = 0, and (K~)* = K+, i.e., 

oo 1 

Kf(x) = T V -[c nPl (Tn) e*™ - s n p 2 (^n) e ±inx ] e~ ne . 
L — ' m 

n=l 
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We will need the commutators of the latter operators. A straightforward 
computation yields 

[K-{x),K+(y)]=C e+ef {x-y) (9) 

where C e {r) = X]^Li( c n el " r + s n e_mr ) e~ ne /n. Inserting Eq. (H) and expand- 
ing 1/(1 — q 2n ) in geometric series we obtain 

C £ (r) = -log[6 e (r)exp (ir/2)] 

with 

. oo 

6 B (r) = -2ie- £ / 2 sin(^^) JJ[1 - 2g 2n e" e cos(r) + q in e~ 2e ] . (10) 

n=l 

Note that 6 e (r) is essentially a regularization of the above-mentioned Jacoby 
theta function, &o( r ) = const. i?i(r/2). We now can define normal ordering, 

x -iVAK e (x) x _ -iv^K+(x) _-h/AKT(x) 

X C X C Kj * 

and a simple computation shows that this amounts to a multiplication with 
the constant b2 £ (0)~ x ^ 2 - This completes our construction of the anyon model. 
One can now compute all anyon correlation functions by using the Hausdorff 
formula and Eq. @, for example the function 

F e ' £ '(x; y) = F £ ' £ \ Xl , . . . , x N ; y h . . . , y N ) : = 

= (n,<f> e (x N y ■ ••(j)e{xi)*^e'{yi) ■ ••</»£' {vn)Q) = 

= Ui< j<k <N b 2e(xk - x j ) x b 2s ,(y j - y k ) x 

which will play an important role further below. 

Having set the stage, we now can describe the second quantization of the 
eCS Hamiltonians if at in Eq. ([[J) and how this leads to a remarkable identity 
which will be the starting point for our solution algorithm: we found a self- 
adjoint operator Ti on T such that the commutator of Tt with a product of 
N anyon operators, 

9?{x) = &M-~ttM, (12) 

is essentially equal to H^^^ f ('x) where is the eCS Hamiltonian with a 
coupling constant determined by the statistics parameter as follows, 

7=2A(A-1). (13) 
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To be more precise: this operator TC obeys the relations 

[w,^(x)]n~^«^(x)n (14) 

with H% as in Eq. ([!]) but with V(r) replaced by the regularized potential 
V2 £ (r) = — d 2 log b 2e (r)/dr 2 , and '~' means 'equal in the limit e [ 0'. We 
note that it is surprisingly simple to construct this operator TC by following 
the arguments for the Sutherland model fll2"] . It has the form 

TC = -—^ ^2 x p(m + n)p(—m)p(—n) * + 

+ (1 - A) J2 n M- n )pi( n ) + p2(-n)p 2 (n)] + ... (15) 

n>0 

where the dots indicate less important terms proportional to QHo, Q 3 and Q 
0. It is interesting to note that the relations in Eq. flip) can be established 
without using Eq. (R), i.e., the second quantization TC obeying Eq. (|TJ) exists 
for a much larger class of Hamiltonians given in Eq. ([!]) with V(r) = Vo(r), 

= E - ( einr + s ' t einr + e_inr ]) e ~ n£ • ( 16 ) 

n=l 

However, as we will see, to obtain an solution algorithm one also needs 

<n,[W,^(x)*^(y)]n> = (17) 

(note that this identity is trivial in the Sutherland case where TCQ = 0, but 
this no longer holds in general): we proved that this identity holds true if 
and only if 

C m C n S m+n ~ S m S n C m+n Vm, 71 — 1,2,... (18) 

which restricts us to c n and s n as in Eq. (||) and thus interaction potentials 
which are Weierstrass elliptic functions (this proof is by a straightforward 
computation using the explicit formula for TC given above ||). We now 
compute the vacuum expectation value of the (trivial) identity 

[Ti, $f (x)*<^ (y)] = -[TC, $f (x)]*d>* (y) + $f (x)*[W, (y)] 

using Eq. (|17j). Inserting Eq. (|I~4D twice we obtain 

#f(x)F £ ' £ '(x;y) * ^'(y)F £ - £ '(x;y) (19) 
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with F e < £ '(x;y) = < ft,$f (x)*$^(y)ft > the anyon correlation function 
defined and computed in Eq. (|TT|) above and regularized cCS Hamiltonians 
acting on different variables as indicated (the bar means complex conjuga- 
tion). This is our remarkable identity and a main result of this letter. 

We now show how Eqs. (|T9|) and ([□]) can be used to construct eigenfunc- 
tions and eigenvalues of the eCS Hamiltonian. The idea is to take the Fourier 
transform of Eq. (|19|), i.e., apply to it (2tt)~ n J d N y e lP ' y (the integration is 
over — 7r < i/j < Ti, of course), and then take the limits e, e' [ 0. To determine 
the possible values for the Fourier modes Pj we observe that b £ (r) x changes 
by a factor e T17rA under r — ► r ± 2ir. Thus the function F £ ' e (x; y) changes 
by a factor e - 17T ( 2N - 2 3+ 1 ) x under yj — > yj + 2tc. The Pj need to be such that 
e lP ' y F e ' e (x; y) is periodic, which enforces Pj = rij + (N —j + |)A with integers 
rij. These are (essentially) the 'quasi-momenta' known from the Sutherland 
model (up to a trivial shift due to a center-of-mass motion contained in our 
eigenfunctions; see Eq. (|33|) belowQ). With that we obtain 

oo 

H N F{x; n) = £ (n)F(x; n) - 7 ^ ^ n 

j<k n=l 



<F(x; n + nE jfc ) + <F(x; n - nE ji; ) (20) 



with (Ejk)t = 8ji — 8k£ for £ = 1, . . . , N, and 

N N 



^o(n)=^P/ = ^[n i + (iV-i + i)A] 2 , (21) 



where the first term on the r.h.s. in Eq. (|20| ) comes from the derivative term 
in the eCS Hamiltonian and partial integration, and the second term comes 
from the interaction terms which we evaluated using Eq. ([DjD. The function 
F(x; n) is the e, e' j 0-limit of the Fourier transform of F e,£ (x; y), i.e., 

F(x; n) = P(x; n) A(x) e^ iAx ^ (22) 

where 

A(x) = Y[b (x k -xj) x (23) 

i<fc 



^^The interested reader can check that removing this center-of-mass motion amounts to 
changing the quasi-momenta as follows, Pj — > n 3 + A(A^ + 1 — 2j) /2, and by relabeling 
rtj — ► riAr+i-j one obtains Sutherland's expression H. 
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and the symmetric, periodic function 



with b e (r) = (1 - e ir " e ) [J^Li [1 - 2q 2n e~ £ cos(r) + g 4n e" 2e ]; the last factor 
in Eq. (p2|) describes an uninteresting center-of-mass motion. We defined 
the functions P as a particular regularization of singular integrals, but these 



integrals can be computed explicitly by expanding in plane waves Hl2 |. Ex- 
panding also in powers on q 2 one obtains V = X]^o ^ ^ wnere the V 1 are 
symmetric polynomials, i.e., their expansion in plane waves has only a finite 
number of non-zero terms (this number of terms diverges as i goes to infinity, 
however) . 

Writing /x = J2j<k^jk^jk with integer fij k and identifying the set of all 
such /x with 1, N< y N ~ i y 2 we now make the following ansatz for an eigenfunction, 

^(x) = 5>(^)F(x;n + M) (25) 

(we suppress the common argument n of ip, a, £ in the following). Then the 
equation 

H N iP = Ei) (26) 

implies 

oo 

[£ (n + /*) - 5] a(/x) = 7 J] n[e£a(/x - nE ifc ) + 4«(/£ + nE ifc )] . (27) 

j<fc n=l 

To solve these equations we make the ansatz 



a(£) = ? M > £ = 5>g^. (28) 

£=0 fcO 

Using s 2 n = T,m=i Q 2nm and c n = 1 + 4 we get 

[£ (n + /£) - £b] «K/£) - ^2 ^m^-m(/£) = 

m=l 

00 

= 7^^7K^(/x-nE jfc ) +7^ S 

J<fc n = l j<k n,m>0 

nm<£ 

xn[«<-„ m (ft - nE jk ) + o^_ nm (/x + nE jfe )] . (29) 
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Eqs. (fZ^)-(P§|) constitute our algorithm to solve the eCS model. We can 
restrict ourselves to (£, /x) such that 

fi jk > -i Vj < k (30) 

[i.e. for other (£, /x) we set a^(A*) = 0] an d determine the c^(^x) and Si from 
Eq. ( |2"9"D recursively. It is important to note that this recursive procedure has 
triangular structure, i.e., amounts to diagonalizing triangular matrices. To 
see that we observe that there is a natural partial ordering, (£', fj,') < (£, /j,) 
if ~ ~ 

£' < £ or (£' = £ and fi' jk < fi jk Vj < k and /x' ^ /x), 

and for fixed (£, /x) there is only finitely many (£',fJ>') with (£',/-i') < (£, /x). 
For n = and £ = we get £o = £o( n ), and «o(Q) remains a free parameter. 
For £ > and /x = we get an equation determining Se as a sum of finitely 
many terms depending only on the a^(/x') and £^ with (£', fj,') < (£, 0), and 
a^(0) remains undetermined. For non-zero /x, there are two different cases. 
If there is a resonance, i.e., if the factor 

£ (n + /x) - £ (n) = 2 [n,- - n fe + (A; - j)A] /i jfc + 

vanishes, we get a linear equation constraining the previously undetermined 
a£>(/j,') for (£',fJ,') < {£,l-i), and a^(/x) remains undetermined. In the generic 
case, i.e. if there is no resonance, a^(/x) is determined as a finite sum of 
terms depending only on ai'(fi') and £# for (£',/-i') < (£, /x). We thus can 
determine all a^(/x) recursively, and at each order £ there remains only one 
of them undetermined [generically this will be a;^(0)]. We note that these 
free parameters correspond to the freedom of choosing different normaliza- 
tions of the eigenfunctions, i.e., changing them amounts to multiplying the 
eigenfunction i/j(x) by a normalization constant which is a power series in q 2 . 
Moreover, it is known that in the Sutherland case one only needs to consider 
TV-tuples n such that 

m > n 2 > • - • > n N > (32) 



since these already provide a complete set of eigenfunctions |L2[ |j, and we 
expect this is true also for q ^ 0. 
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It is clear that resonances make our algorithm somewhat more involved, 
and it is therefore interesting to mention some cases where resonances can 
be ruled out. For example, there is never a resonance for fi > and n as in 



Eq. (|32|) ||12|| , and therefore resonances do not occur in the Sutherland case 
q = 0. Moreover, for N — 2, it is easy to see that resonances can only occur if 
A is integer. However, for iV > 2, there are infinitely many resonances which 
are independent of A, e.g. for N = 3 and n such that rii — 2n 2 + n 3 = 3u 
with integer u, one has a resonances for all fi such that /X13 = — v — {i\ 2 and 
fi 2 3 = 2v + (/ii2 arbitrary integer), and for rational values of A, additional 
'coincidental' resonances (i.e., they depend on A) are to be expected. Thus, 
for N = 3, resonances can be ruled out for irrational A and n such that 
(rii — 2n 2 + n ? ^)/3 is non-integer. Obviously, a more general analysis of the 
occurrence and implications of resonances would be welcome. 
The eigenfunctions which we get are of the form 

</>(x) = J n (x|g) A(x) e^ iiVA ^ (33) 

where we defined J n (x|g) = YlT=o ^n( x ) Q 2e with 

i 

J n(x) = ^^a^(M)^'(x;n + ^). 

£'=0 fj. 

It is interesting to note that the latter sums are always finite, i.e., one can 
prove highest weight relations for the functions V which imply that there 
are only finitely many \i obeying Eq. flSTf ) and such that V (x; n + n), < 
£' < £, are different from zero |K|. Moreover, as already mentioned, the 



J n (x|g) are uniquely determined up to normalization. In the case q = our 



algorithm reduces to the one in Ref. [12|, and the results there imply that the 



J n (x|g = 0) = Jn( x )' n obeying the condition in Eq. (|3"2l), are proportional 
to the Jack polynomials [|]. It is therefore natural to regard the J n (x|g) as 
elliptic generalization of the Jack polynomials. 

Recently it was shown that the formal power series expansion of the 
eigenfunctions of the eCS model in q 2 converge in the L 2 -norm sense pi 



This result suggests that the formal power series which we obtained actu- 
ally converge, and in particular, that our elliptic generalizations of the Jack 
polynomials are well-defined symmetric functions. For q = it is known that 
they define a complete orthogonal set of eigenfunctions 0, and we conjecture 
the same is true also for finite q. Obviously, a more detailed investigation of 
these functions would be welcome. 
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We conclude with a two remarks. Firstly, we note that the only result 
from our anyons construction which we actually needed to get our solution 
algorithm is the remarkable identity in Eq. ([19]). For e = e' = the latter 
reduces to the following identity of elliptic functions, 



N 



d 

E(^-^) F ( x -y) = Y. 2X ^- l A p{Xi ~ Xk) ~ p{y ^ yk) ] F ^ ]y) (34) 

j=i X i y i j<k 
where 

„, , Ui< j<k < N MU x k- x j}) x MUyj-yk\y 

F(X| Y) = ' llS^MIfe^ ■ <35) 

We suggest that it is this remarkable identity which distinguishes the quantum 
Calogero-Moser-Sutherland models [rather than the existence of a product 
ground state function of the form as in Eq. (fffi) ] since it can be used to 
obtain an algorithm to solve the model and holds for the elliptic well. 
After this work was completed we found an elementary, direct proof of this 
based on the following identity of Weierstrass elliptic functions, 

[({x) + C{y) + C(z)} 2 = p{x) + p{y) + p{z) iix + y + z = U 

(whose proof is given as an exercise on p. 446 in Ref. ||13|| ) HTo| . Still, the 
generalization which we got from the anyons contains important additional 
information: the e-regularization is needed to fix a unique self-adjoint exten- 
sion of the differential operator in Eq. ([I]). Finally, we note that it would 
be interesting to find an argument to obtain from the identity in Eq. (|i~9|) an 
algorithm providing finite series representations of the eigenfunctions, e.g., 
by finding a way to solve Eq. (p7|) non-perturbatively (i.e., without expanding 
in q 2 ). 



Acknowledgements: 

I thank A. Carey and A. Polychronakos for their interest and helpful 
discussions and K. Takemura for pointing out to me Ref. ||21|| . A helpful 
remark of K. Fredenhagen is acknowledged. This work was supported by the 
Swedish Natural Science Research Council (NFR). 



11 



References 

[1] Calogero, F., J. Math. Phys. 10 (1969) 2191 and 2197; and 12 (1971) 
419. 

[2] Sutherland, B., Phys. Rev. A4 (1971) 2019 and A5 (1972) 1372. 

[3] Moser, J., Adv. Math. 16 (1975) 197; Calogero, F., Lett. Nuovo. Cim. 
13 (1975) 411. 

[4] For review see Olshanetsky, M.A. and Perelomov, A.M., Phys. Rept. 71 
(1981) 313 and 94 (1983) 313. 

[5] See e.g. Polychronakos, A. P.: "Generalized statistics in one dimension," 
|hep-th/9902157| . 

[6] See e.g. Jain, S., Mod. Phys. Lett. All (1996) 1201; Guhr, T., Muller- 
Groeling, A., and Weidenmuller, H.A., Phys. Rept. 299 (1998) 189. 

[7] Forrester, P.J., Nucl. Phys. B388 (1992) 671 and J. Stat. Phys. 72 
(1993) 39. 

[8] Macdonald, I.G.: Symmetric functions and Hall polynomials. Oxford 
Mathematical Monographs. Clarendon Press, Oxford 1979; Stanley, 
R.P, Adv. in Math. 77 (1989) 76. 

[9] Langmann, E.: 'Second quantization of the elliptic Calogero-Sutherland 
model', math-ph/0 10205. 

[10] Langmann, E.: 'An algorithm to solve the elliptic Calogero-Sutherland 
model', (to appear on math-ph). 

[11] Awata, H., Matsuo, Y., Odake, S., Shiraishi, J., Phys. Lett. B347 

(1995) 49 and Nucl. Phys. B449 (1995) 347; Iso, S, Nucl. Phys. B443 
[FS] (1995) 581; Marotta, V., and Sciarrino, A., Nucl. Phys. B476 

(1996) 351. 

[12] Carey, A.L., and Langmann, E., Commun. Math. Phys. 201 (1999) 1. 

[13] Whitaker, E.T., and Watson, G.N.: Course of modern analysis, J^th 
edition, Cambridge Univ. Press, 1958. 



12 



Etingof, P.I., and Kirillov, A.A., Duke Math. J. 74 (1994) 585; Ruijse- 
naars, S. N. M., J. Math. Phys. 40 (1999) 1595. 

Felder, G., and Varchenko, A.: "Integral representation of solutions 
of the elliptic Knizhnik-Zamolodchikov-Bernard equations," Int. Math. 



Res. Notices 5 (1995) 221 Qhep-th/9502165P and "Three formulas for 



eigenfunctions of integrable Schroedinger operators," |hep-th/951112ti| . 

see e.g. Azuma, H., and Iso, S., Phys. Lett. B331 (1994) 107; Iso, S., 
and Rey, S.J., Phys. Lett. B352 (1995) 111; Yu, Y., Zheng, W., and 
Zhu, Z., Phys. Rev. B56 (1997) 13279. 

Wen, X.G., Phys. Rev. B41 (1990) 12838 and Phys. Rev. Lett. 64 
(1990) 2206. 

Gradshteyn, I.S., and Ryzhik, I.M.: Table of integrals, series, and prod- 
ucts, 5th edition, Academic Press, Boston, 1994. 

Reed, M., and Simon, B.: Methods of modern mathematical physics. II. 
Fourier analysis, self-adjointness, Academic Press, New York, 1975, p. 
177. 

Carey, A.L., and Hannabuss, K.C., J. Func. Anal. 75 (1987) 128. 

Takemura, K, Lett. Math. Phys. 53 (2000) 181; Komori, Y., and 
Takemura, K.: "The perturbation of the quantum Calogero-Moser- 
Sutherland system and related results," |nath . qa/0009244 . 



13 



